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We show that the process of spontaneous symmetry breaking can trap a field theoretic system in a 
highly non-trivial state containing a lattice of domain walls. In one large compact space dimension, 
a lattice is inevitably formed. In two dimensions, the probability of lattice formation depends on the 
ratio of sizes L x , L y of the spatial dimensions. We find that a lattice can form even if R = L y /L x 
is of order unity. We numerically determine the number of walls in the lattice as a function of L x 
and L y . 
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I. INTRODUCTION 

As a system undergoes spontaneous symmetry break- 
ing, it relaxes into a new phase. If the topology of the 
vacuum manifold is non-trivial, there can be topological 
defects that slow down the relaxation process. Eventu- 
ally, however, it is believed that the system will get to its 
ground state. It is then somewhat surprising to find that 
certain systems will essentially never get to the ground 
state. Instead they will be trapped in a highly non-trivial 
state which, in one compact spatial dimension, consists 
of a lattice of domain walls. 

In the present paper, we build upon a large body of 
earlier work. Analysis of how simple domain wall net- 
works relax may be found in Refs. 0, 0, 0, 0- The kind 
of systems we are interested in lead to more complex do- 
main walls and are discussed in Refs. 0, IE IS El EU 
ITU IT^ . Il3| . Most of the discussion so far has been in 
one spatial dimension. Here we study a phase transition 
in two dimensions in which a network of domain walls is 
produced. Then we study the evolution of the network 
toward its final lattice state. We focus on the likelihood 
that the system does not reach the vacuum, and on the 
properties of the final lattice state. 



II. THE MODEL 



We start with the SU(5) x Z 2 model: 
L = Tr(d M $) 2 - V{$) 



(1) 



where $ is an adjoint scalar field represented by 5 x 5 
Hermitian, traceless matrices. The potential is assumed 
to satisfy V($) = V(— $) and the Z 2 symmetry is under 
$ — > — If we truncate the model to just the four di- 
agonal matrices, the only surviving symmetry is the per- 
mutation of the diagonal elements of $. Now the model 



contains four real scalar fields (i = 1,..,4), and the 
Lagrangian is : 
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where the fields fi are defined by: 

= h(x)X 3 + f 2 (x)X s + f 3 (x)T 3 + U(x)Y , (4) 

where A3, As, r 3 and Y are the diagonal generators of 
SU{5): 
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= -diag(l,-l,0,0,0) 
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=diag(l, 1,-2, 0,0), 
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:diag(2,2,2,-3, -3) 



(5) 



The model now has an 5 5 x Z 2 symmetry corresponding 
to permutations of the 5 diagonal elements of the SU (5) 
and the Z 2 transformation corresponding to change of 
sign. The potential is minimized by a vacuum expec- 
tation value of fi and the symmetry is broken down to 
S 3 xS 2 . 



The symmetry breaking: 

S 5 x Z 2 -» S 3 x S 2 



(6) 



leads to (5!/3!2!) x 2 = 20 distinct degenerate vacua. 
Each vacua is labeled by the vacuum expectation value 
of $ and these will be denoted by: 



±(2,2,2,-3,-3) 
±(2,-3,2,2,-3) 
±(2,2,-3,-3,2) 
±(-3,2,2,-3,2) 
±(-3,2,-3,2,2) 



±(2,2,-3,2,-3), 
±(-3,2,2,2,-3), 
±(2,-3,2,-3,2), 
±(2,-3,-3,2,2), 
±(-3,-3,2,2,2) 



We can find a domain wall solution interpolating be- 
tween any two distinct vacua. Hence there are a large 
number of domain walls in the model. But there are two 
distinct sets of domain walls (Fig : the first set con- 
sists of domain walls across which the vacua differ by 
the Z-i transformation, while the second consists of do- 
main walls that separate vacua that do not differ by Z^. 
The 10 different Z% walls can further be distinguished 
according to their masses Q. There are 3 solutions that 
are least massive. As an example, one of these light- 
est walls interpolates between vacua in the directions 
(2, 2, 2, —3, —3) and —(—3, —3, 2, 2, 2) i.e. vacua related 
by 2 permutations and a sign. Then there are 6 interme- 
diate mass walls between vacua related by 1 permutation 
and a sign e.g. ±(2,2,2,-3,-3) and -(-3,2,2,2,-3), 
and 1 heaviest wall interpolating between vacua related 
by permutations and a sign e.g. +(2, 2, 2, —3, —3) and 
-(2,2,2,-3,-3). 

The charge of a wall (up to normalization) is defined by 
the difference $(+oo) — $(— oo). For the lightest Z 2 walls, 
the charges are: ±(— 1, — 1, 4, — 1, — 1) and permutations 
of the entries. Hence even the lightest Zi wall comes in 
5 different charges (position of the 4 entry) and there are 
5 corresponding antiwalls. We shall say that a wall and 
antiwall are of the "same type" if their charges only differ 
by a minus sign. Otherwise the walls are of "different 
type". 

In Ref. [llj it was shown that walls and antiwalls of 
different type repel in this model. This led to the possi- 
bility of a stable lattice of alternating domain walls and 
antiwalls provided that all the neighboring walls and an- 
tiwalls are of different type. If a neighboring wall and 
antiwall are of the same type, they will attract, then an- 
nihilate. The lattice can be schematically depicted as in 
Fig. |21 or by the sequence of vacuum configurations: 



+ (2, 2, 2, -3, -3) -» -(2, -3, -3, 2, 2) 
+(-3,2,2,-3,2) -» -(2,-3,2,2,-3) 
+ (2, 2, -3, -3, 2) ^-(-3, -3, 2, 2, 2) 
+ (2,2,2,-3,-3)^... 



(7) 



This sequence gives the minimal lattice of 6 kinks. 

In our one dimensional simulations [l3T | we found that 
a domain wall lattice is inevitably formed after a phase 



S 5 /S 3 xS 2 




FIG. 1: The vacuum manifold of the model consists of 10 + 10 
points. Domain walls correspond to solutions interpolating 
between any two vacua. There are two kinds of domain walls, 
one interpolating between the two sets of 10 vacua related 
by Z2 transformations, and the other interpolating between 
vacua within a set of 10 vacua. If we fix the vacuum at x = 
—00 there are 10 ways of choosing the vacuum at x — +00 
that differs by a Z2 transformation, leading to 10 different Z2 
walls. Of these 10 Z2 walls, 3 are of lowest energy_j_6 are of 
intermediate energy, and 1 has the highest energy 
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FIG. 2: Depiction of the sequence of vacuum configurations 
in the domain wall lattice. 



transition provided the compact dimension is not too 
small. In our preliminary work in two equal compact spa- 
tial dimensions (L y — L x ) we had found that a domain 
wall network with junctions is formed after the phase 
transition and this network coarsens, eventually ending 
in the vacuum. In other words, a domain wall lattice was 
not seen to form. The present paper deals with the case of 
two unequal compact spatial dimensions (L y 7^ L x ). We 
expect that when one of the dimensions is very small, the 
system will behave like the one dimensional case, and a 
lattice will almost certainly be formed. When the sizes 
of the two dimensions are comparable, based on the pre- 
liminary results, we would not expect a lattice to form. 
This will turn out to be incorrect and our more extensive 
analysis will show that even in the L y — L x case there is 
a small chance of lattice formation. We will be interested 
in determining the domain wall lattice characteristics as 
a function of the relative sizes of the dimensions of the 
compact space. 
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III. NUMERICAL TECHNIQUE 

In our simulations we try to mimic some features of a 
typical non-equilibrium second order phase transition. In 
order to do that we use a Langevin-type equation based 
on the fundamental Lagrangian of the theory Eq. @: 

(df-V^fi + diV + Ddtf^Ti, (8) 

The Langevin equation describes the evolution of the sys- 
tem coupled to a heat bath with temperature T. The ef- 
fect of the bath is modeled by including in the equations 
of motion a dissipation term with dissipation constant D, 
and a stochastic force Ti(x,t). The stochastic force is a 
Gaussian distributed field obeying: 

(r i (a; ) t))=0, 

{T, t {x,t)T 3 {x' ,t')) ^2DTS lJ S{x- x')S(t - t') (9) 

The amplitude of the noise in Eq. J^J is chosen so as to 
guarantee that independently of the initial field configu- 
ration and of the particular value of the dissipation, the 
system will always equilibrate toward a thermal distribu- 
tion with temperature T. 

In a typical simulation, we start with an arbitrary ini- 
tial condition and evolve Eq. © for a period of time long 
enough for the system to equilibrate. We then mimic an 
instantaneous phase transition to T = temperature, by 
setting the noise term in Eq. © to zero. In the subse- 
quent period of evolution, the dissipation forces the fields 
to its possible vacua values, and a complex network of 
domain walls forms, separating different vacuum regions. 
This network is then allowed to evolve until a final sta- 
ble state is reached. A detailed discussion on how the 
evolution proceeds can be found in [l3j . 

The equations of motion were discretized using a stan- 
dard leapfrog method and periodic boundary conditions 
were used. The model parameters were set to m = 
1/(2V6), A = 1/2 and h = -3/40. For this particular 
choice of parameters there is a known explicit analytical 
solution for the relevant domain profile [3, l| ■ The lattice 
spacing was Sx = 1. and the time-step set to 5t = .5. For 
the parameters above, the wall core is resolved by more 
that 10 lattice points which should be accurate enough 
for the desired purposes. The value of the dissipation 
coefficient D does not influence the results during the 
stochastic stage of the simulation and we set it to D = .5 
to ensure rapid thermalization. We keep the same dissi- 
pation value during the evolution regime at T = 0, the 
reasons behind this choice and its effects on the final re- 
sults are discussed in detail in the next Section. 

In order to identify the walls during and at the end 
of the simulation, we convert the fi(x)'s back into their 
original SU(5) x Z 2 matrix form, $(ar), via Eq. Q. It is 
easy to see that for the lowest energy kinks in the model, 
the sign of Tr[$ 3 ] vanishes at the defect core ^J. For 
these kinks, one (and only one) of the fields <&jj crosses 
zero at the core as well. In the simulation we measure, at 
each time-step, the value of Tr[$ 3 ] and look for pairs of 



neighboring lattice points for which this quantity changes 
sign. In these cases we determine which of the field com- 
ponents goes through zero. This way we are able both to 
find the domain walls and to identify them according to 
their five possible charges, as described in Section ITT1 

All the simulations were run in two-dimensional do- 
mains with periodic boundary conditions. For each set 
of runs, the size of the simulation box in the x-direction 
was fixed to a constant L x , whereas its y dimension, L y 
was allowed to vary. For each choice of L x and L y , we 
executed several independent runs with different random 
seeds and averaged the relevant physical quantities over 
those to obtain the final results. 



IV. RESULTS 

In Fig. |3 we have an example of how the evolution of a 
typical network of walls proceeds in a toroidal domain. In 
this case, the larger dimension was set to L x — 500 and 
the smaller to L y = 150. For early times after the tem- 
perature is quenched to zero (top plot, t = 500), we can 
observe a dense population of walls forming a complex 
network. While it is not obvious from the topology of the 
vacuum manifold that domain wall junctions ("nodes") 
exist in the model 0, the simulation clearly shows their 
presence. These nodes may be viewed as being due to 
the existence of domain wall lattices in one dimension. 
To see this, imagine a circle in two spatial dimensions. It 
is possible to have a lattice of walls on this circle. Now 
if we shrink the circle to a point, the walls forming the 
lattice have to converge to a point. In other words, the 
walls on one of the circles are like spokes on a wheel and 
have to necessarily have a point of convergence in two 
dimensions. While an infinite number of lattice configu- 
rations can be constructed, the minimal number of walls 
in a stable lattice is six. In the simulation also we see that 
a majority of intersecting nodes consist of six incoming 
walls. As the evolution proceeds, several mechanisms 
play a role in decreasing the overall wall density. Walls 
and anti-walls of the same type annihilate each other and 
closed loops can be formed later decaying into radiation. 
The intersecting nodes are also dynamical and can dis- 
entangle when nodes and anti-nodes collide. A snapshot 
of the wall network for intermediate times can be seen 
in the middle plot in Fig. |3| Finally, for long times, all 
the loops collapse and all nodes annihilate, leaving be- 
hind a set of straight parallel walls wrapped around the 
shortest dimension of the simulation torus. Depending 
on their charge, some of these walls may still annihilate 
each-other, and the final state of the evolution is either 
the vacuum or a lattice of equidistant repelling walls. In 
the bottom plot of Fig. [3] we can see that for t — 5000 the 
system has relaxed into a stack of parallel walls. Since 
we know how to identify the wall charge from the simu- 
lation data, once this stage of the evolution is reached it 
is straightforward to predict which walls will annihilate 
and how many (if any) will remain in the final lattice. In 
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50 100 150 200 250 300 350 400 450 500 

FIG. 3: The domain wall network for three typical stages in an 
evolution in a toroidal domain, with dimensions L x — 500 and 
L y — 150. The different shades correspond to the 5 possible 
charges the domain walls can have in the theory. Note that 
in the bottom figure there is a pair of neighboring wall and 
anti-wall of the same type (the walls just before and after the 
300 mark). These will later on annihilate each other leading 
to a final stable lattice consisting of 10 walls. 
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FIG. 4: Probability of forming a lattice of domain walls as a 
function of the ratio of the two torus dimensions. L x is fixed 
for each curve, taking the values 500 (circles), 200 (squares) 
and 100 (triangles). Standard deviation error bars (not shown 
for clarity) are of the same order of magnitude as the result. 



this particular case there is only one wall anti-wall pair 
that annihilates, resulting in a lattice of 10 walls. 

Though the pattern of the evolution remains the same 
independently of the relative dimensions of the two- 
dimensional torus, its final outcome depends markedly 
on the values of L x and L y . 

The most striking trend in the data is that the final 
state of the evolution is very likely to be a stable wall 
lattice when L y <C L Xl whereas the vacuum state is the 
predominant outcome when the two dimensions are com- 
parable. This dependence is illustrated more explicitly in 
Fig. 21 where we show the probability of having a lattice 
forming as a function of the ratio of the two dimensions 
of the simulation domain. For each choice of L x and 
L y , the formation probability is defined as the number of 
independent runs that led to the formation of a lattice, 
divided by the total number of runs. For the three fixed 
choices of the largest dimension, L x = 100, 300, and 500, 
the results were obtained by averaging over 64, 40 and 40 
runs for each value of the smallest dimension L y . In all 
three cases, L y was varied from to a larger value, close 
to L x in the two cases of small L x but only to about 250 
for L x — 500 (the reason being that for the larger case 
it takes a very long time for the system to reach its final 
state, and simulations for higher L y become numerically 
demanding). The errors calculated from the standard 
deviation of the data are of the same magnitude as the 
results. We chose not to show them in the plot (and 
likewise in the following figures) so not to obscure the 
results. 

For L y = 0, the probability of forming a lattice is quite 



high, being in fact equal to unity in the case of L x = 
500. This is to be expected since this case corresponds 
to having the evolution takin g pl ace in a one-dimensional 
domain. As was shown in 11], in ID, as the size of 
the physical domain goes to infinity, the probability of 
forming a lattice tends to one. On the other hand, for 
the symmetric 2D system where L x — L y , the lattice 
formation probability is quite small. In the case of the 
L x = L y — 100 simulations, out of 64 independent runs 
only 3 had a stable lattice of walls as their final state. 
A similar result was found for the L x = L y — 200 case, 
where 3 lattices were formed in a total of 32 simulations. 
We also observed that in the few cases where a lattice is 
formed in a L x = L y box, the lattice shows no preferred 
direction of alignment 1 . In the low and medium range 
of L y / L x on the contrary, the walls in the final lattice 
are always parallel to the y direction, as expected (see 
Fig. |3J). For L x — 200 the first lattices in the x direction 
appear only for L y = 150 whereas for L x = 100 these 
are observed for values of the y dimension larger than 
L y = 70. 

It is interesting to note that the formation probability 
for large values of the ratio R = L y /L x seems to be 
independent of the choice of L x . A similar effect can be 
seen in the averaged number of walls in the final state, 
as plotted in Fig. |3J The curves for the three choices 



In the L x = L y case, the domain wall lattice spontaneously 

breaks the spatial symmetry under interchange of x and y in 

every member of the ensemble but the choice of direction by 
which it is broken is random. 
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FIG. 5: Average number of walls in the final stable configura- 
tion. The spatial dimensions of the torus are the same as in 
Fig. 01 Standard deviation error bars (not shown for clarity) 
are of the same order of magnitude as the result. 



FIG. 6: Total wall length per unit area in the final stable 
configuration. The spatial dimensions of the torus are the 
same as in Fig. 0] Standard deviation error bars (not shown 
for clarity) are of the same order of magnitude as the result. 



of L x converge for a value of roughly L y /L x = 0.4 and 
coincide above that. Other quantities seem to follow the 
same pattern in this regime. The size ratio above which 
walls parallel to the larger dimension are found in the 
final state, is given by R = 150/200 = 0.75 in the case 
L x = 200. For the smaller box L x = 100, the result 
is similar with walls in the x direction appearing above 
R = 70/100 = 0.70 

The outcome of the evolution clearly seems to vary 
within two well defined regimes. For R > 0.4 the prop- 
erties of the final state depend only on the ratio of the 
dimensions of the compact space, whereas for lower val- 
ues of R the results depend also on L x . In the small R 
regime, for example, the final number of walls varies con- 
siderably. That this is to be expected can be easily seen 
by considering the L y = limit, corresponding to the 
one-dimensional case. In this situation we would expect 
the final number of walls to be proportional to the length 
of the ID domain. In other words, the final wall density 
p - defined as the total length in walls divided by the 
area of the lattice - should depend little on L Xl which is 
indeed observed: 

Pico = 0.047, p 200 = 0.046, p 5 oo = 0.042 (10) 

As shown in Fig. the curves for the number density 
are quite close to each other for small values of L y /L x , 
diverging for higher ratios of the two dimensions. 

Finally we note that, as mentioned in Section IIIII a 
non-zero dissipation term is kept throughout the whole 
evolution. We expect that the presence of dissipation, 
though changing the type of scaling of the dynamics of 
the network, should not influence qualitatively the out- 
come of the evolution. The reason for this is that the 
probability of forming a lattice should depend mostly on 



the number of uncorrelated domains in the initial state of 
the evolution. As a check, we ran a few simulations with 
dissipation set to zero during the evolution stage. The fi- 
nal results turned out to be different from the above, but 
only for small values of L x and R. The reason behind 
this is that after the quench there is still a lot of energy 
stored in the field which can be re-distributed among its 
different degrees of freedom. One consequence of this 
process is that it can effectively destroy the initial pat- 
tern of vacuum domains. For example, in the L x = 500, 
L y = case, the formation probability is still unity but 
the average wall number in the final state is less than 
shown in Fig. 01 This observation can be interpreted as 
a consequence of the effective domain size being changed 
soon after the quench, when the dissipation is turned off. 
This would also explain why the effects vanish for larger 
boxes. The conclusion is that keeping the dissipation on 
during the evolution allows us to use smaller simulation 
boxes to probe the large scale limit of theory. 

V. SUMMARY 

We have numerically studied phase transitions in a 
model with S5 x Z2 symmetry breaking down to S3 x S2 
in two compact spatial dimensions. As expected, the 
outcome depends on the ratio of the sizes of the two di- 
mensions: R = L y /L x . As R — > 0, the system behaves 
as in one spatial dimension and a domain wall lattice 
forms. As R — ► 1, the system behaves as in two spatial 
dimensions and a domain wall lattice rarely forms. What 
is surprising though is that the critical value of R below 
which a lattice forms with large probability is of order 
unity: R c ~ 0.4. 
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In earlier work [13( we have studied the evolution of 
the network in time and found that the total length of 
the network falls off as t~ 0J in the L y = L x (R = 1) 
case. Here we have focused on the late time state of the 
system with various values of R and have determined the 
domain wall density p at late times as a function of R. At 
the moment we do not know if there are naturally occur- 
ring systems in which a domain wall lattice exists. The 
domain wall lattice may also be relevant for cosmology, 
possibly in the context of brane cosmology. 
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